Free Energy Barrier for Electric Field Driven Polymer Entry into Nanoscale Channels 
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Free energy barrier for entry of a charged polymer into a nanoscale channel by a driving electric 
field is studied theoretically and using molecular dynamics simulations. Dependence of the barrier 
height on the polymer length, the driving field strength, and the channel entrance geometry is 
investigated. Squeezing effect of the electric field on the polymer before its entry to the channel is 
taken into account. It is shown that lateral confinement of the polymer prior to its entry changes 
the polymer length dependence of the barrier height noticeably. Our theory and simulation results 
are in good agreement and reasonably describe related experimental data. 

PACS numbers: 87.15.A-, 36.20.Ey, 87.15. H- 



Intcraction of polymer chains with nanoscale chan- 
nels is a phenomenon of rich basic science and numer- 
ous potential applications [H-0] . Better understanding of 
related biological processes, analyzing biopolymers like 
DNA Q , investigating the existing theories on static and 
dynamic properties of confined polymers [3jL and sep- 
aration of polymers of different lengths [iHH are sam- 
ple advantages of the experiments in which polymers are 
forced to interact with nanochannels. Also, this phe- 
nomenon provides possibility of observation and direct 
study of protein-nucleic acid interactions Q and nucleic 
acids secondary structures A polymer entering into 
a narrow channel experiences an entropic barrier due to 
the reduction of its conformations and an external force 
should be applied to overcome this barrier. A voltage 
difference in the case of charged polymers like DNA [8|, 
or flow injection may provide the driving force 0, [13] • 

The entry process of a charged polymer into a 
nanochannel in a wall under a voltage difference can be 
divided into two stages 



(a) 





First, the polymer comes 



to the channel vicinity by a pure or electrically biased 
diffusion and one of its ends finds the channel entrance. 
Then, if the driving field overcomes the entropic bar- 
rier, the polymer enters the channel. Calculation of the 
barrier height as a function of the system parameters 
such as the driving field strength, the channel geome- 
try, and properties of the polymer is of great importance 
for numerous related experiments. Polymer translocation 
through nanochannels and polymer escape from entropic 
traps which recently has been used as a polymer sepa- 
ration method 0, [1| are examples of such experiments. 
Despite extensive experimental and theoretical studies 
of driven polymer threading into nanoscale channels and 
numerous suggestions for enhancing the polymer capture 



rate [ 1 lM 1 3l . |17| , the entry barrier height and its depen- 



dence on the system variables is less investigated. 

In this paper, we calculate free energy barrier expe- 
rienced by a flexible charged polymer entering into a 
nanochannel under an applied electric field theoretically 
and using coarse grained molecular dynamics (MD) sim- 



FIG. 1: Schematic of a polymer under the electric field of 
strength E behind a channel of diameter d in a wall, (a) 
The polymer can be viewed as a 2D chain of blobs of size Lj. 
(b) A polymer which is also confined laterally by a cylinder 
of diameter D. The polymer can be considered as a chain of 
blobs closely packed inside a cylinder of height Ln. (c) A 
segment of the polymer entered to the channel inside which 
the electric field is of strength Ei. 



ulations. We obtain dependence of the barrier height 
on the polymer length and the driving field strength for 
two different geometries of the channel entrance. The 
external driving field could have squeezing effect on the 
polymer before its entry to the channel and affect its 
static and dynamic behaviors. Our study shows that this 
squeezing effect and lateral confinement of the polymer 
prior to the channel entry noticeably affect the depen- 
dence of the barrier height on the polymer length and 
the field strength. Our theoretical and MD simulation 
results support each other and are in reasonable agree- 
ment with some related experimental data. 

To obtain free energy barrier for two channel entrance 
geometries shown in Fig. [TJ we calculate the polymer free 
energy before its entry to the channel and after entering 
a fraction of its monomers. We assume that electric fields 
inside and outside of the channel are uniform and their 
strength may be different. As case /, consider the poly- 
mer which contains N charged monomers of diameter b 
and charge q, pushed to a wall by the electric field of 
strength E Q (see Fig. [I] (a)). We assume that there is no 
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N qEb/k T ' NqEb/k T 



FIG. 2: (color online) (a) /3 in L oc and 7 in .R||L-4 oc TV 7 versus the polymer length. With increasing TV crossover from 
values corresponding to case I to those of case II can be seen. (b),(c) Simulation results for polymer entry barrier and fitted 
functions of the theory for cases I and II. Coefficients a,— 3.02, b=1.50 and c=0.32 are fitting parameters. Considering scaled 
axis of these figures, agreement between simulation and theory can be seen. Insets: entry barriers versus E" 1 and fitted lines. 



electrostatic interaction between monomers (correspond- 
ing to high salt concentration) and that the monomers 
only interact with the electric field. The polymer de- 
creases its electric energy by spreading on the wall in 
price of loosing a part of its entropy. The polymer layer 
thickness, Li, can be determined by balancing entropic 
and electric energies. In this layer, the polymer can be 
viewed as a 2D chain of blobs of size Lj, where the num- 

ber of monomers in each blob is gx ~ (nf) " and v is the 
Flory exponent. Entropic energy of the polymer is k B T 

per blob, F ent ~ k B TN Cj~jj ■ F° r case H-i consider 
a similar polymer which is also confined from sides by a 
cylinder of diameter D (see Fig. Q] (b) ) . Because of simul- 
taneous confining effects of the electric field and the cylin- 
der, this is similar to a polymer in a closed space which 

its blobs are closely packed. The number of monomers in 

i 

each blob is given by gu ~ (-^jt) 3 "~ 1 , where f2 ~ LjjD 2 
is the confinement volume. The entropic energy of the 

polymer in this case is F ent ~ k B T ( ^rr~ ) 1 @- 

Electric energy of the polymer layer, F e i c , in both of 
these cases with its zero taken on the wall is ~ NqE L m 
(m = /, II) . The equilibrium layer thickness can be ob- 



tained by minimizing total free energy F 
with respect to L m : 



F, 



,fqE b\ 1 +" h fqE Q b\ (Nb 2 ^ 3 " 



k B T I \D 2 



Number of monomers per blob in the two cases are 
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(1) 



(2) 



Free energy per monomer behind the channel is ~ ksZ. . 

Consider a polymer in each of the two cases that is 
going to enter a channel of diameter d on the wall. 
Entropic energy of each monomer after entry into the 



channel is ~ in which, gi n = (f) 1 ^ is the num- 

ber of monomers per blob inside the channel. Elec- 
tric energy change of n monomers after entering to the 
channel is ~ —nqEil, in which Ei is the electric field 
strength inside the channel and I = -^-d = n(^) 1 /"d 
(see Fig. [2(c)). Accordingly, total free energy change 
of the polymer after entry of n monomers is AF m ~ 

k B T -±)n- rPqEifyV'd. Free energy bar- 

rier height for the polymer entry, AU m , is the maximum 
of AF m with respect to n and is given by 
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qE t d 



(3) 



In simulations for checking our theory, we model 
the polymer by N spherical monomers of diameter a 
connected by FENE potential. Shifted and truncated 
Lennard- Jones potential is used to model monomer- 
monomer and monomer-wall excluded volume interac- 
tions and it is assumed that E Q = Ei = E. All our simu- 
lations are performed with ESPResSo 14] using Langevin 
thermostat to keep the temperature fixed. 

To check our theoretical results on static behavior of 
the polymer before its entry to the channel, we measure 
some variables in our simulations by keeping the chan- 
nel closed and letting the polymer to fluctuate behind it. 
Lateral gyration radius of a polymer confined between 

two planes of spacing L is Ru ~ Nib^^ 4 [9J. For a 
fixed diameter of the confining cylinder, D = 12a, ex- 
ponents /? and 7 in relations L oc N 13 and R\\Lt- oc N 1 
are obtained via local slopes of the log-log plots of L and 
versus N. Results for two different field strengths 
are shown in Fig. [2] (a). For a fixed value of D, short 
and long polymers correspond to cases / and 77, respec- 
tively. With increasing N, crossover from case I with 
expected exponents /?/ = and 77 = 0.75 to case II 
with (3[i = ~ 0.6 and -fu — ~ 0.14 can be seen 
in Fig. [5] (a). The value of a in relation L oc E~ a is 
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AF JV 
k B T ~ gn 

experimental result of Ref. 



also obtained from simulations for three different poly- 
mer lengths. For both cases I and 77, it is found that 
a ~ 0.4, in agreement with the theory. Results on static 
behavior of the polymer is of importance in the subject of 
polymer compression by external force in nanochannels 
15]. As an example and for comparison with experi- 
ment, the change in the free energy of a polymer arising 
from its confining in a narrow channel under an elec- 
tric field can be written from case II of our theory as 
oc A 1 ' 6 . This is in excellent agreement with 

To check the main result of our theory, Eq. [3j we mea- 
sure the barrier height for the polymer entry for each 
given set of N and E in both of cases I and II. We hx 
the end monomer of the polymer at the channel entry 
and leave the other monomers to equilibrate under the 
applied electric field. Then we release the polymer end 
to probe whether its entry is successful or no. We obtain 
~~^T as ^ ne logarithm of successful entry probability 
(fraction of at least 10 3 simulations in which the poly- 
mer entry is successful). Results for polymers of different 
lengths at different values of the electric field strength are 
shown in Fig. |3](b) and (c). As it can be seen, there is a 
very good agreement between our theory and simulation 
results in both of the two cases. In the insets of Fig. 
[5] (b) and (c), deviation of AU versus E~ x from linear 
form in large intervals of E in both cases, no dependence 
of entry barrier on polymer length in case I contrary to 
case II, and decrease of entry barrier with the polymer 
length in case II can be seen. These results are different 
from predictions of the previous theories: linear decrease 
of the barrier height with the applied voltage [13| and its 
dependence on the polymer length even for simple pores 



17| . In previous studies, dependence of the barrier height 
on N is attributed to the restriction of the polymer end 
to the channel during its entrance [l8j |. Because of the 
compression effect of the driving field which is taken into 
account in our work, only monomers in the first blob feel 
the restriction and its energy cost is independent of N. 

One should note that the validity range of our theory 
in both cases is 1 < g m < N. g m > N means that 
the electric field is too weak to has confining effect on 
the polymer behind the channel and cannot overcome 
the entropic barrier. Consistent with this description, 
in simulations with smaller polymers and lower electric 
fields, we obtained a ~ and /3 ~ v. g m < 1 means that 
the polymer loses most of its entropy behind the chan- 
nel and feels no entropic barrier on its entry (as simula- 
tion results also show in Fig. 2 (b) and (c)). Existence 
of two voltage limits, the one below which no polymer 
entry happens and another one above which voltage de- 
pendence of the capture time changes from exponential 
to linear function [ly, Qj| could be related to these limits. 

In voltage driven polymer translocation experiments, 
one of the measurable quantities is the capture time, r, 
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FIG. 3: (color online) Logarithm of capture time versus AV 
extracted from referenced articles (symbols) and fitted equa- 
tion of our theory, Eq. [4] (solid lines). All the data sets are 
from polymer translocation experiments through a-hemolysin 
channel. As it can be seen, all data sets are followed well by 
the suggested function. Fit parameters are shown in Tabled] 



TABLE I: Parameters £, 77, and A obtained from fitting of Eq. 
13] to seven data sets shown in Fig. [3] 





Experimental data 


c 


V 


A 


1 


ssDNA, cis side (Ref. [12]) 


-5.42 


3.14 


0.60 


2 


ssDNA, cis side (Ref. [20]) 


-4.33 


2.38 





3 


ssDNA, trans side (Ref. [20]) 


0.27 


3.22 


1.77 


4 


ssDNA, cis side (Ref. [21]) 


-3.01 


3.11 


0.90 


5 


DS, cis side (Ref. [19]) 


-7.09 


2.35 


0.01 


6 


DS, trans side (Ref. [22]) 


-6.13 


3.85 


0.51 


7 


PSS, cis side (Ref. [23]) 


-6.01 


2.83 






at different voltage differences, AV" [H E3, • The 

average capture time can be written from Van't Hoff- 

Arrhenius law as r = tq exp (j^pf^j [20], in which AU 

is the free energy barrier height, t^ 1 is the frequency 
of polymer fluctuations behind the channel when entry 
of folded polymer is possible 24{ or frequency of the 
polymer ends reaching to the channel entrance when en- 
try of folded polymer is not possible. From Van't Hoff- 
Arrhenius law and eq. [3] of the theory, and considering 
that l/(l + i/) ~ ~ 0.6, the equation 



lnr = C+ ry - A 



( qEjb 
\k B T 



qE t b 



(4) 



can be used for fitting to experimental data by 77, and 
A as the fit parameters, regardless of the similarity of the 
channel entrance geometry to case I or II of the theory. 
In equation |4] we have assumed Ei ~ =j- and E a = 5Ei, 
where £ is the channel length and S < 1 (there is no satis- 
factory knowledge about values and spatial functionality 
of Ei and E a ). Also, b and q are obtained by course- 
graining the polymer such that the size of effective spher- 
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ical monomers is equal to the polymer diameter. For ex- 
ample, for translocation of single stranded DNA (ssDNA) 
through a-hemolysin channel at the room temperature, 
we use t = lOnm, b — Iran, and q = 3e (e is the elec- 
tron charge), which results in ~ O.OlAV(mV). It 
should be noted that from the Van't Hoff-Arrhenius law 
C = In r and from Eqs. [2] and El T) ~ (f)*"* 7 for 

1 | 1 

both of the two cases, A ~ (|) 2 2u 8 0m6 for case I and 
A ^ (^+^o.6 (i^- 6 for case jj 

In Fig. [3l lnr versus Ay reported for translocation 
of three flexible polymers, ssDNA, dextran sulfate (DS), 
and poly(styrenesulfonic acid) (PSS), from cis and trans 
sides of a-hemolysin channel and fitted function of Eq. 
[4] are shown. The values of fitting parameters are shown 
in Table HI C m Eq. |4] depends only on To which is out 
of interest here, rj depends on the ratio of the channel 
diameter to the monomer size and as all the seven data 
sets shown in Fig. [3] are for experiments with the same 
channel, a-hemolysin, similar values of r\ is reasonable. 
The most important parameter which is a measure of 
the squeezing effect of the electric field on the polymer 
before its entry to the channel is A. For data sets 2, 5 
and 7, where polymer is added to the cis side and applied 
voltages are relatively low (up to 160 mV), A is negligible. 
This is because of the weak electric fields which results in 
small squeezing effects. Nonzero value of A with similar 
applied voltages for the trans side (data sets 3 and 6), 
could be because of the smaller diameter of the channel 
in this side which causes the electric field in the vicinity 
of the channel entrance to be stronger and the parameter 
6 to be larger [51]. It should be noted that according to 
ref. [l2j ]. compression of the polymer inside the cis side 
vestibule of a-hemolysin channel happens very rarely. So 
entry of the polymer from both the cis and trans sides 
is similar to case I of the theory (and not the case II). 
For the cis side, data sets 1 and 4, applied voltages are 
higher (up to 280 mV) which makes the squeezing effect 
more considerable and the value of A nonzero. 

In summary, free energy barrier height for electric field 
driven polymer entry into channels of two different en- 
trance geometries has been calculated theoretically and 
by MD simulations. The barrier height dependence on 
the electric field strength and the polymer length, the 
effect of lateral confinement of the polymer before its en- 
try to the channel on these dependencies and relevance 
of results with experiment have been discussed. 

We are grateful to Mohammad R. Kolahchi for useful 
comments. 
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